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Abstract. The luminous high-redshift (z — 2.17) quasar 
S5 0836+710 has been observed at 5 GHz with the VSOP. 
We compare the properties of three images obtained 
from the observation: a low-resolution ground array im- 
age (dynamic range 4600:1), a full- resolution VSOP im- 
age (900:1), and an image made with only the space base- 
lines (200:1). The space baselines alone are sufficient for 
a reliable recovery of the source structure, within the lim- 
its of the achieved spatial sampling of the visibility data. 
The curved jet ridge line observed in the images can be 
described by Kelvin-Helmholtz instabilities developing in 
a relativistic outflow with the Mach number of about 6. 
This description holds on the scales of < 700/i -1 pc, and 
is shown to be consistent with variable apparent speeds 
observed in the jet. 
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1. Introduction 

The VSOP (VLBlQ Space Observatory Program) is 
a Space VLBI (SVLBI) mission utilizing the world- 
wide array of radio telescopes and an orbiting 8-meter 
antenna deployed on the Japanese satellite HALCA0 
( Hirabayashi 1996). The satellite has an elliptical orbit, 
with the apogee at 21000 km, perigee at 560 km, and or- 
biting period of about 6 hours. In each observation, the 
data stream from the satellite is recorded by a network 
of 5 STS (Satellite Tracking Station), and subsequently 
correlated with the data from participating ground tele- 
scopes. Regular VSOP observations started in Septem- 
ber 1997 at 1.6 and 5 GHz. A 5 GHz VSOP observation 
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of 0836+710 (4C71.07, z=2.17) was made on October 7, 
1997. This radio source is an ultraluminous quasar with 
well established, correlated (Otterbein et aL 1998) broad- 
band variability in gamma-ray (Fichtcl ct al. 1994), X- 
ray (Brunncr ct al. 1994), optical (von Lindc ct al. 1993), 
mm- and cm- radio regimes (Marscher & Bloom 1994). 
The source has a well collimated VLBI-scale jet (Hummel 
et al. 1992) extending out to ~180 milliarcseconds (mas). 
The jet is substantially curved, with lateral displace- 
ments of its ridge line and oscillations of its transverse 
width (possibly correlated with the observed jet speeds, 
Krichbaum ct al. 1990). The observed speeds change re- 
near the core, apparent speeds 
= 0.5 and H n = 



markably along the jet: 
Aipp ~ 10 hr 1 c are measured (with qo 



100/ikms Mpc ); at 3 mas core separation, the speed 



decreases down to /3 a 



2-3 hr 



and then becomes 



larger once again, further out (Krichbaum et al. 1990, Ot- 
terbein 1996). This kinematic behavior may be caused by 
MHD instabilities developing in the jet plasma (Hummel 
et al. 1992). 

We present here a 0.2-mas (0.8/i _1 pc) resolution 
VSOP image of 0836+710, and discuss the morphology 
and properties of the compact jet in this quasar. 

2. VSOP observation and data reduction 

0836+710 was observed with the VSOP at 5 GHz for 11.5 
hours, with the VLBA^ array (see Zensus, Diamond, & 
Napier 1995, and references therein) providing ground 
support for the observation. The data were recorded in the 
VLBA format, using total observing bandwidth of 32 MHz 
divided in two intermediate frequency (IF) bands, each 
having 256 spectral channels. The STS in Usuda (Japan), 
Tidbinbilla (Australia), and Robledo (Spain) were used 
for the satellite data acquisition. The data were corre- 
lated at the VLBA correlator in Socorro, with output pre- 
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0836+710 at 4.962 GHz 1997 Oct 07 
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Fig. 1. The u?>coverage the VSOP observation of 
0836+710. The inner circles correspond to the baselines 
between ground radio telescopes; the space baselines are 
grouped into larger ellipses extending at a P.A.w 45°. 
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Fig. 2. Visibility amplitudes (top panel) and phases 
(bottom panel) on the baseline between HALCA and 
VLBA:North Liberty, after amplitude calibration, fringe 
fitting, and phase self-calibration. The lines drawn 
through the data represent the final CLEAN model of 
0836+710. 



averaging time of 1.966 and 0.524 seconds for the ground 
and space baselines respectively. Fringe visibilities were 
detected in the HALCA data recorded at Tidbinbilla and 
Roblcdo. The resulting sampling function (wy-coverage) of 
the final correlated dataset is shown in Figure [l], indicating 
an improvement of resolution by a factor of ^3 compared 
to ground VLBI observations at the same frequency. 

2.1. Fringe fitting 

Post-processing of the correlated data has been done in 
AIPS0 and DIFMAP ( [Shepherd et al. 1994|) . We applied 
amplitude calibration, using the antenna gain factors and 
system temperature measurements. After inspecting the 
IF bandpasses, the last 46 channels were flagged in each 
IF, owing to significant (50-80%) amplitude reduction. 
This has reduced the total observing bandwidth to 26.2 
MHz. We applied the phase-cal information available for 
the VLBA antennas, and then corrected the residual de- 
lays and rates, using the single-band (SB) and multi-band 
(MB) delay fringe fitting. We used solution intervals of 2 
(SB) and 3 (MB) minutes, and accepted all solutions with 
SNR > 7. After the fringe fitting, the residual phase vari- 
ations were found to be within 3° on the ground baselines, 
and < 10° on the space baselines. We then averaged over 
all frequency channels, and calibrated the phases with a 
point source model (to enable time averaging). Finally, 
the data were exported into DIFMAP and further time- 
averaged into 60-second bins. The amplitude and phase 
errors were calculated from the scatter in the unaver- 
aged data. Figure || gives an example of typical HALCA 
baseline visibility data. The estimated RMS noise on the 
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HALCA baselines is about 4.5 times higher than on the 
ground baselines. 

2.2. Imaging 

To investigate the impact of the VSOP mission, we have 
imaged both the entire dataset (VSOP dataset hereafter) 
and only the ground baseline data (VLBA dataset), and 
compared the properties of these images. Both images 
have been produced after several cycles of hybrid map- 
ping and self-calibration. 

We use natural weighting for gridding the ground array 
data, enhancing the sensitivity to extended emission at the 
price of slightly decreased image resolution. The gridding 
weights are also scaled by amplitude errors raised to the 
power of —1. Both phase and amplitude self-calibration 
have been used, with amplitudes being allowed to vary 
only after the total model flux has approached the zero 
spacing flux to within 3%. The image obtained using this 
procedure is shown in the left panel of Figure 

To image the VSOP dataset, we use uniform weighting 
which provides a better angular resolution at the expense 
of lowering slightly the sensitivity to extended structures. 
Because the estimated noise on HALCA baselines is much 
higher than on the ground baselines, scaling the gridding 
weights by the amplitude errors weights down significantly 
all of the long iw-spacings. To avoid this, we switch off 
the amplitude scaling of gridding weights. As a safeguard 
measure, we apply only phase self-calibration to the VSOP 
dataset. After a good fit to the data is achieved, we ad- 
just the antenna gains, correcting for small constant off- 
sets between the model and the visibility amplitudes. An 
example of the obtained fit to the data can be found in 
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Fig. 3. Ground array (left) and VSOP (right) images of 0836+710, with the jet ridge line marked. In the ground 
array image, the dotted line rectangle shows the area covered by the VSOP image. The image parameters are as 
follows. The ground array image: restoring beam is 2. 15x1. 75 mas at P.A.= 30.8°; contour levels are drawn at —0.7, 
0.7x 1.5" mJy/beam (n = 0,...,17); the peak flux density is 918mJy/beam. The VSOP image: restoring beam is 
0.93x0.31 mas at P.A.= —36.3°; contour levels are drawn at —1.5, 1.5x1.5™ mJy /beam (n = 0, 13); the peak flux 
density is 432mJy/beam. 



Figure The final high-resolution image is shown in the 
right panel of Figure [| 

To test the reliability of the satellite data, we make an 
additional image using only the space baselines (HALCA 
image) . The resulting image shown in Figure is produced 
with only phase self-calibration applied. The structures 
seen in Figure || are consistent with the structures found in 
both VSOP and VLBA images in Figure |[ This reassures 
the good quality of the data received from the orbiting 
antenna. 

Table compares the properties of all three images. 
Note the roughly 5 times smaller dynamic range of the 
VSOP image, which is caused by the poorer sampling and 
increased noise of the data on the space baselines. The suf- 
ficiently high dynamic range of the HALCA image (200:1) 
reflects once again the high quality of the satellite data. 



Table 1. Parameters of the images 



Image 


Stot 


5max 




CRMS 


T-jpcak 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


VLBA 


2.167 


0.918 


-0.001 


0.2 


4600:1 


VSOP 


2.156 


0.432 


-0.003 


0.5 


900:1 


HALCA 


0.731 


0.325 


-0.007 


1.4 


200:1 



Column designation: 2 - total CLEAN flux [Jy]; 3 - peak flux 
density [Jy/beamj; 4 - lowest negative flux density [Jy/beam]; 
5 - RMS noise [mJy/beam]; 6 - peak-to-RMS dynamic range. 

3. Properties of the compact jet 

In both VSOP and VLBA images, the jet is continuous, 
with several enhanced emission regions embedded in it. In 
Figure || curvature of the jet is noticeable even at < 1 mas 
scales, and it is consistent with the images obtained at 
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Fig. 4. Image of 08364-710 obtained using only the 
baselines to HALCA. All structures seen in the image 
are consistent with those found in both the VSOP and 
VLBA images. This shows the reliability of the data 
recorded on the space baselines. The restoring beam is 
0.68x0. 18 mas at P.A.= —36.9°; contour levels are drawn 
at -3.2, 3.2xl.5"mJy/beam (n = 0, ...,11); the peak flux 
density is 325mJy/beam. 
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Fig. 5. Relation between the the ridge line and observed 
jet speeds in 08364-710. Circles are measured offsets of the 
jet ridge line from the jet axis at P.A.= —146°; squares 
denote speeds reported by Otterbein (1996) at different 
locations in the jet; solid line represents apparent speed 
variations derived from the measured ridge line offsets, for 
a jet with Mach number, Mj = 6, Lorentz factor, 7j = 11, 
and opening half-angle, <fij = 1°. 



higher frequencies, most notably at 86 GHz (Otterbein et 
al. 1998). We consider now the observed curvature of the 
jet represented, in the maps in Figure by the curved 
ridge line of the jet. We use the median P.A.= —146° 
to approximate the direction of the jet axis, and mea- 
sure the absolute values of corresponding offsets of the 
jet ridge line. The offsets are plotted in Figure || (opaque 



circles). One can see that the offsets show a remarkable 
periodicity, possibly with several periods superimposed. 
The derived offsets also appear to be correlated with the 
apparent speeds measured along the jet (Otterbein 1996, 
black squares in Figure |5|) . 

The apparent curvature of the jet may reflect the pres- 
ence of plasma instabilities in the jet (e.g. MHD or Kelvin- 
Helmholtz instabilities). The amplitude of ridge line dis- 
placements from the jet axis can be used to estimate the 
jet Mach number, Mj, using the analytical approximation 
obtained by Hardee (1984) for the amplitude, Q(r), of the 
fastest spatially growing mode of Kelvin-Helmholtz insta- 
bility: Q(r) = G(ro)(r/ro) e . Here r refers to the distance 
along the jet axis, "0" denotes the initial values, and e 
is determined by the jet Mach number and opening half- 
angle, (f>j, so that e = 0.83/ (Mj sin^j). 

To quantify the initial conditions of the jet plasma, we 
adopt the results from Otterbein et al. (1998) obtained 
by modelling spectral evolution in the core of the VLBI 
jet in 08364-710. The initial amplitude of instability can 
then be approximated by the size of the core, so that 
G{ r o) = 0.8pc, at r — 15pc corresponding to the es- 
timated distance from the jet nozzle to the core. With 
the reported jet bulk Lorentz factor, 7j = 11, viewing 
angle, 6j = 3.2°, and opening half-angle </>j = 1°, we 
obtain the Mach number Mj ss 6, within the innermost 
3 mas of the jet. It is similar to the value reported in 
ptterbein 1996 , and we take this as an argument in favor 
of our choice of the direction of the jet axis. We expect the 
inner jet to be pressure confined, with external pressure 
P ext » 25Pj e t- Farther down the jet, at distances > 10 mas 
(> 700 /i _1 pc), the derived values of Mj appear to increase 
up to ~ 30 — 50, if we assume that the jet viewing angle 
remains the same. 

Perhaps a more attractive possibility is to preserve the 
Mach number along the jet, and instead use variations of 
the jet viewing angle for explaining the observed ridge line 
offsets. For the same set of Mj , 7j , and <j>- } , we determine the 
variations of the viewing angle required to reproduce the 
offsets. The derived 6 } changes smoothly between 3° and 
65°, with largest values occurring at ~ 3 mas and ~ 10 mas 
distances — which are incidentally the locations where 
the most pronounced lateral displacements are reported 
(Krichbaum et al. 1990, Otterbein 1996). The combina- 



tion of the assumed 7j and derived 8j allows us to predict 
the apparent speed, /3 a pp, along the jet. The predicted /3 a pp 
(plotted in Figure [| with the solid line) is in a reasonable 
agreement with the measured jet speeds. At r > 13 mas, 
the predicted /3 ap p becomes flat, as it approaches the limit 
of (7? — l) 1 / 2 . At distances r > 15 mas, a moderate in- 
crease of 7j up to ~ 15, could provide a better agreement 
with the measured speed. However, at these distances, the 
measurement uncertainties for both the speeds and ridge 
line offsets are too large to warrant arguing in favor of 
variable bulk speeds along the jet. We would like to em- 
phasize that the achieved consistency between the model 
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and observed speeds is preserved, if the chosen direction 
of the jet axis are within a few degrees from the one we 
have used. In that case, only the derived jet Mach number 
would have differed from Afj = 6 obtained above. Sparsity 
of the available speed measurements precludes better kine- 
matic modelling and selecting the jet axis unambiguously. 
We expect that making such a selection would require de- 
termining the continuous velocity distribution along the 
jet, for which a suitable monitoring program should be 
designed, aimed at obtaining a sequence of high dynamic 
range and high fidelity VLBI images of the source. 

4. Conclusions 

The reported observation of 0836+710 shows that the 
VSOP mission provides an excellent opportunity for high- 
resolution and high dynamic range imaging of VLBI-scale 
jets in extragalactic objects. The obtained VSOP image 
of 0836+710 allows to investigate the jet morphology and 
physical conditions on linear scales up to 1 kpc. We have 
studied the curved ridge line in the jet, and shown that, 
on the scales of <700/i _1 pc, it can be modelled by a 
pressure confined relativistic outflow with Lorentz factor 
7j « 11, Mach number Mj w 6, and opening half- angle 
<pj w 1°. In this description, the two main factors respon- 
sible for the observed curvature of the jet ridge line are: 
1) Kelvin-Helmholtz instabilities developing in the rela- 
tivisitic plasma and 2) variations of the angle, 0j, between 
the velocity vector of the outflow and the line of sight. We 
have shown that the derived variations of 0j can be used 
for predicting the jet kinematics, and are consistent with 
the observed variable apparent speeds along the jet. 
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